Transition to linear domain walls in nano-constrictions 
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Domain walls in nano-constrictions are investigated with a focus on thermal properties. In general, 
the magnetization component perpendicular to the easy axis which in a domain wall usually occurs 
has a value different from the easy-axis bulk magnetization value with a separate phase transition 
at a critical temperature below the Curie temperature. Since this effect is the more pronounced 
the smaller the domain wall width is we investigate it especially in domain walls with a confined 
geometry, using analytical arguments, mean-field theory and Monte Carlo simulations. Our findings 
may contribute to the understanding of magneto-resistive effects in domain walls with sizes of only 
a few atomic layers, as e. g. in nano-contacts or nano-constrictions. 
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Well controlled domain walls could become important 
constituents of future magneto-electronic devices |l| . Es- 
pecially, the understanding of domain walls in confined 
nanometric geometries is important since those can show 
behavior deviating from their usual bulk properties, as 
e. g. a controlled pinning and a strongly reduced do- 
main wall width 0, HH The latter is thought to 
contribute to large magneto-resistance effects of domain 
walls in nano- wires [^, nano-constrictions 0, and nano- 
contacts 

In only few publications the temperature dependence 
of domain wall properties was investigated [10, ■ 
Most important in this context is the pioneering work 
of Bulaevski and Ginzburg ^3 '^ho showed within the 
framework of Ginzburg-Landau-theory that for a one di- 
mensional domain wall profile (e. g. a Bloch wall) the 
easy-axis and hard-axis components of the magnetiza- 
tion, respectively, are two separate order parameters with 
different critical temperatures. In other words, the per- 
pendicular magnetization component which arises neces- 
sarily in a domain wall has at finite temperatures val- 
ues lower than the easy-axis equilibrium magnetization 
(leading to the term "elliptical domain walls") and van- 
ishes completely for a temperature Th which is lower 
than the Curie temperature of the bulk material (lead- 
ing to the term "linear domain walls" for temperatures 
Th<T<T,). 

However, the deviation of from Tc is proportional 
to the squared inverse domain wall width 13] and, hence, 
should be very small. Consequently, linear walls are hard 
to detect experimentally [lj|, their relevance is consid- 
ered to be rather low and most of the numerical cal- 
culations of domain wall properties are performed using 
micromagnetic codes where the assumption of a constant 
magnetization value is made in contradiction to the find- 
ings described above. In the following we will investigate 
in how far thermodynamic deviations from pure Bloch- 
like domain wall structures can become relevant due to 
the reduced size of domain walls in confined geometries. 
This is an important question since it was suggested 
that linear domain walls might explain the observed large 



FIG. 1: Sketch of the wall geometry. The magnetization of 
the first and the last plane are antiparallel along the easy axis. 
The magnetization rotates in a plane via the hard axis. 



magneto-resistance effects in constrained domain walls. 

In the following we consider a domain wall structure 
as shown in Fig. ^ In the picture, each arrow represents 
the mean magnetization of a plane. Our Monte Carlo 
(MC) simulations are for a full three dimensional model 
allowing for magnetization fluctuations within the planes 
while in the mean-field (MF) approximation each (infi- 
nite) plane is represented as one magnetization vector so 
that the model is effectively one dimensional. Depend- 
ing on the details of the methods we apply we use either 
fixed boundary conditions where the first and the last 
plane of the system are fixed as shown in the figure or we 
use anti-periodic boundary conditions. In both cases we 
force a domain wall into the system which for large sys- 
tem size (number of planes N) will not fill up the whole 
system. However, for smaller system size the boundary 
conditions force the domain wall to adopt the system size. 
This models a domain wall caught in a nano-constriction 
2] where perfect pinning is assumed. 

We investigate the system in terms of a classical spin 
model with spin variables \Si\ — 1 on a cubic lattice 
with lattice constant a and energy contributions from 
ferromagnetic exchange between nearest neighbors with 
coupling constant J and a uniaxial anisotropy with D > 
defining the easy axis of the system. 
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In the zero temperature limit all spins within a plane 
will be parallel and the domain wall profiles can be calcu- 
lated analytically in the continuum limit where the en- 
ergy density (per cross-sectional area) is a one dimen- 
sional integral 
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The magnetization of the domain wall rotates within a 
plane and can be expressed by a single angle of rotation 
9. The Euler-Lagrange-equation which minimizes the en- 
ergy above is solved by an elliptic integral 
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where Sq — a\J J /2D is the zero-temperature domain 
wall width of an unconstrained wall and the integration 
constant c is given by the boundary condition x{t:/2) — 
L/2 with L = {N — l)a. The wall profiles are then given 
by Jacobian sine and cosine functions 



sn 



/ex 1 



Sh{x) 



cn 



j'CX 1 



(2) 



The limit c —> 1 corresponds to an unconstrained wall 
with the usual Bloch-wall profiles 

Se{x) — tanh(x/(5o) Sh{x) = cosh~^{x/So). (3) 

The opposite case, c » 1, corresponds to a very con- 
strained wall {L So) where the wall is forced to adapt 
the system size and the profiles follow simple trigonomet- 
ric functions, 

Se{x) = sin(7r2;/L) Sh{x) = cos(7ra;/L). (4) 

In this limit the actual domain wall width is L. 

As an example. Fig. [3 shows a comparison of these 
expressions numerical MF results, obtained as described 
below in the low temperature limit /cbT = 0.02J. The 
smaller system shows the case of a very constrained wall 
where the wall profiles are already described by trigono- 
metric functions while the larger system shows an inter- 
mediate case. Note that in Ref. corresponding cal- 
culations were made, but for other boundary conditions 
and system geometries, respectively. 

However, we want to focus on thermal properties and 
in order to obtain results for finite temperatures we start 
with the MF Hamilton operator which (neglecting terms 
without S) is 



Hmf = 
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Hm,+m,+i) + D{S^f (5) 
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FIG. 2: Ground state easy- and hard-axis magnetization pro- 
files for constrained walls. Data points are from low tem- 
perature MF calculations, solid lines correspond to Eq. |5| 
D/J = 0.003. 
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FIG. 3; Easy- and hard-axis magnetization profiles in MF 
approximation for two diS'erent temperatures. D/J = 0.03 



where the contribution ASi ■ rrii comes from the four 
neighbors within the plane. Then we solve the MF self- 
consistency equations. 



m,: = (S,) = iTrS,;e' 
-'HyiF/kBT 
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with Z — Tve ' numerically. Here, mi is the 

thermally averaged magnetization of the ith plane and 
the trace is an integral over the unit sphere. These equa- 
tions can be solved iteratively, starting with an arbitrary 
magnetization profile and then let the equations evolve 
until a stationary state is reached. 

Let us start with the case of an unconstrained wall. 
Fig. 131 shows domain wall profiles for two different tem- 
peratures where the system size L = 40a is large enough 
so that an equilibrium domain wall for an anisotropy 
value oi D = 0.03J fits well into the system. The solid 
lines are the analytical functions as calculated above, 
here, simply the tanh and 1/cosh profiles. Interestingly, 
for finite temperatures the mathematical form of the wall 
profile is conserved, solely the amplitudes and the domain 
wall widths vary with temperature. Hence, in the limit 
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FIG. 4: Temperature dependence of easy- and hard-axis order 
parameters and (reduced) inverse domain wall width from MF 
calculations. L = 50a, D/J = 0.3. 
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FIG. 5: Anisotropy dependence of the two critical tempera- 
tures Tc and Th- Data points are from MF calculations as in 
Fig. 21 solid hues from Eq. El 



i S> (5o the thermodynamic wall profiles can be described 
as in Eqs. El but with a temperature dependent domain 
wall width 5{T) and temperature dependent amplitudes 
Mh{T) and Me(T). Furthermore, Fig. 01 suggests that 
that the temperature dependence of the two amplitudes 
is not the same. 

Instead, these amplitudes, Mh and Mg, define two dis- 
tinct order parameters as shown in Fig. 01 Obviously, 
these two order parameters vanish continously at two 
different temperatures where the upper one is the usual 
Curie temperature Tc and the lower one is a second criti- 
cal temperature, T?i, which describes the phase transition 
of the hard-axis component of the magnetization vector. 
For temperatures Th < T < Tc the domain wall is linear 
with an easy-axis component of the magnetization only. 
For T < Tfi in general one finds an elliptical wall pro- 
file which for lower temperature goes over to the usual 
Bloch wall profile. Also shown in Fig. ^is the reduced, 
inverse domain wall width demonstrating that 6(T) in- 
creases slightly with temperature for elliptical domain 
walls, shows a kink at Th, and then diverges in the linear 
wall regime approaching Tc- 



Both critical temperatures depend on the strength 
of the anisotropy D as is shown in Fig. [51 Tc cor- 
responds to the usual bulk Curie temperature and its 
anisotropy dependence expanded with respect to D/J is 
fceTc/J = 2 + j^D/J + 0[D/JY. As Fig. demon- 
strates, Th has also a linear dependence on D/J in the 
range of anisotropics which is shown in the figure and our 
numerical results suggest k^Th/J « 2 — 0.53D/J. This 
means that the difference between both critical tempera- 
tures scales with the squared inverse domain wall width 



kB{Tc~Th)/J -QA0a^/5l. 
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These findings, the second phase transition, its depen- 
dence on 5q and the diverging domain wall width are 
qualitatively in agreement with the earlier calculation 
within the framework of the Ginzburg-Landau theory [l^ 
— an expansion close to Tc — where it was also shown 
that close to Tc the (linear) wall has a tanh-profile. How- 
ever, for experimental systems reasonable anisotropics 
are rather small, so that in general the two critical tem- 
peratures should nearly coincide and Mh should be close 
to the easy-axis magnetization Mg. Nevertheless, experi- 
mental investigations of linear domain walls exist. In [l^ 
the influence of the wall structure (either elliptical or lin- 
ear) on the domain wall mobility was investigated. Here, 
the deviation of Th from Tc was 1%. However, since — 
as Eq. 1^1 shows — these effects increase with decreasing 
domain wall width larger effects should occur in smaller 
domain walls as found, e. g. in constrained geometries. 

For a strongly constrained domain wall and in MF 
approximation Th can be estimated analytically. For a 
system with antiperiodic boundary conditions and small 
anisotropy, so that L <C the wall profiles follow 
trigonometric functions which means that the angle of 
rotation from plane to plane changes in each plane by 
the same amount -Ka/L. This leads to a decrease of the 
MF coming from the two adjacent planes proportional to 
cos(7ra/L). Including this in the MF Hamiltonian (Eq. 
jsj the critical temperature for zero anisotropy can be 
calculated in the usual way, now leading to 



4 2 ,7ra, 
kY{Th/J = 3 + 3C0s(-^). 



(7) 



Note, that for a constrained domain wall with L <ti Sq 
much bigger effects may occur than before (Eq. EJ. 

In Fig. El we compare the formula above with numeri- 
cal calculations for the case of constrained walls. We use 
D = 0, so that the condition L ^ i5o is always fulfilled, 
and simulate systems with fixed boundary conditions cal- 
culating Th from the hard-axis magnetization. Fig. |H| 
demonstrates that pronounced effects can be found when 
the domain wall is constrained to only a few atomic lay- 
ers. The agreement of our MF data with Eq. |3 is very 
good. Nevertheless, the MF approximation strongly un- 
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FIG. 6; Dependence of the critical temperature Th from the 
size of the constrained domain waU. Comparison of MC data, 
MF calculations and Eq. Q (solid line) . The dashed line is a 
guide to the eye. 



derestimates this size dependent effect: this can be con- 
cluded from the fact that even in the extreme case of a 
three layer system — which due to the fixed boundary 
condition corresponds to a free monolayer — a finite Th 
is found even though it is clear from the Mermin- Wagner 
theorem that no order should occur. 

Therefore, we additionally used MC methods to inves- 
tigate the break-down of ferromagnetic order in the origi- 
nal spin-model (Eq. ^ more rigorously. Using MC meth- 
ods thermal excitations (thermally excited spin-waves) 
are fully taken into account. We use a heat-bath algo- 
rithm and single-spin- flip dynamics for the simulations. 
At every MC step each spin is subject to a trial step 
consisting of a small deviation from the original direc- 
tion The lateral dimension of the system is up to 
128 X 128 with periodic boundary conditions. The num- 
ber of planes is varied from 3 to 21 where we use fixed 
boundary conditions for the first and last planes. We 
start the simulations with an abrupt domain wall and 
let it relax for lOOOOMCS. Then we calculate the abso- 
lute value of the magnetization component perpendicu- 
lar to the easy axis, averaged over the whole system and 
for another lOOOOOMCS as order parameter of the phase 
transition. Note that the precise definition of the order 
parameter is important: without calculating the absolute 
value of the perpendicular magnetization the long time 
average of one magnetization component will always be 
zero since the wall magnetization might rotate in the hard 
plane. Furthermore, to average over the whole system is 
also important since the wall can move diffusively. 

Th is then determined from a finite-size scaling analysis 
of the order parameter above where the lateral system 
size is varied from 8 x 8 to 128 x 128. The resulting data 
points are also shown in Fig. |H| The scaling analysis 
works well with the exponents (3 and v from the three 
dimensional Heisenberg model. Only for a very small 
number of planes deviations occur indicating a crossover 
from three to two dimensional behavior |l7l | . As expected 



Th now goes to zero in the limiting case of a tri-layer 
system so that the change of the critical temperature is 
more dramatic. E. g., for a wall consisting of 5 atomic 
planes {L = 4a) Th is reduced by about 35% as compared 
to the Curie temperature and even for a temperature of 
0.5Tc the hard-axis magnetization (the degree of order 
within the wall) is reduced by about 20% as compared 
to the easy-axis bulk value. 

To summarize, investigating the influence of thermal 
activation on the properties of domain walls in nano- 
constrictions we have demonstrated that with increasing 
temperature Bloch wall profiles change via elliptical walls 
to linear domain walls. The temperature range where 
these effects occur scales with the squared inverse domain 
wall width so that in general it is rather small. However, 
since in confined geometries the relevant quantity is the 
size of the constriction larger effects can be found. The 
break-down of ferromagnetic order is due to the fact that 
in a domain wall the mean exchange field decreases due 
to the finite angle of rotation between neighboring mag- 
netic moments. Hence, it is a general effect which will 
also occur in other types of domain walls as, e. g., vortex 
walls. Our findings may have an impact on the under- 
standin g of domain wall magneto-resistance properties 
[isl ll^Tigl Eoj , especially its temperature dependence, 
for two reasons: first, in an elliptical or linear domain 
wall the degree of spin disorder is larger in the wall than 
in the bulk of the domain since the value of the order pa- 
rameter is lower. Second, as suggested in in a linear 
domain wall the change of the magnetization direction 
is abrupt while only its value is changing. This means 
that the spin of a conductance electron passing a linear 
wall cannot follow a continuously rotating magnetization 
direction. 
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